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Two-loop SUSY QCD correction to the gluino pole mass 

Youichi Yamada a * 

a Department of Physics, Tohoku University, Sendai 980-8578, Japan 

We calculate the pole mass of the gluino as a function of the running parameters in the lagrangian, to 0(o? s ) 
in SUSY QCD. The correction shifts the pole mass from the running mass by typically 1-2 %. This shift can be 
larger than the expected accuracy of the mass determination at future colliders, and should be taken into account 
for precision studies of the SUSY breaking parameters. The effects of other corrections are breifiy commented. 
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1. Introduction 



Extention of the standard model by supersym- 
metry (SUSY), with the breaking scale not much 
higher than the electroweak scale, has been stud- 
ied as a very promising solution to the hierarchy 
problem between the electroweak scale and the 
Planck/grand unification scale. In these mod- 
els, such as the minimal supersymmetric standard 
model (MSSM) [Q, all particles in the standard 
model have their superpartners with masses be- 
low O(l) TeV. These new particles are then ex- 
pected to be produced at colliders in near future, 
such as the CERN Large Hadron Collider (LHC) 
and the International Linear Collider (ILC). 

One of the main motivations for the experimen- 
tal study of these new particles 0, the SUSY 
particles, is the determination of the soft SUSY 
breaking parameters |3I4| . which gives an impor- 
tant information of the SUSY breaking mecha- 
nism in the unified theory. For example, the uni- 
fication of three gaugino masses (M 3 , M 2 , Mi) at 
the same scale as that of the gauge couplings is 
a crucial test for the SUSY grand unified theory 
|5I6| and superstring phenomenology [7|- 

For this purpose, in addition to the precise 
measurements of the physical parameters of the 
SUSY particles 0, one also needs precise pre- 
diction of the relations between these observables 
and parameters in the lagrangian. In some cases, 
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we have to calculate the relations beyond the one- 
loop order to match expected experimental pre- 
cision. For example, the two-loop mass correc- 
tions have been calculated for the top and bottom 
quarks [8], squarks in the first two generations [5], 
and Higgs bosons [TU|. 

Here we focus our attention to the mass of the 
gluino g, the SU(3) gaugino, in the framework of 
the MSSM. At the LHC, gluino is, if it is suffi- 
ciently light, expected to be copiously produced 
Its mass nig can be determined from the 
distributions of the decay products. A study [T5] 
shows that, for the SUSY parameter set SPSla 
given in Ref. [T3| with nig ~ 600 GeV, nig can 
be determined to accuracy dnig = 8 GeV from 
precision data at the LHC with 300 fb -1 , and 
even to Snig = 6.5 GeV when combined with the 
data from the ILC. On the other hand, the one- 
loop QCD contribution to the the difference be- 
tween the pole mass nig and the running mass M3 
of the gluino [14115116] is much larger, typically 
O(10) %. One therefore naively expect the two- 
loop mass correction might be 0(1) %, similar 
to the experimental uncertainty. It is therefore 
important to examine whether higher-order cor- 
rections to the gluino mass is really relevant in the 
determination of the SUSY breaking parameters 
at future precision measurements. 

In this talk we present the pole mass of the 
gluino as a function of the lagrangian parameters, 
including 0{a 2 s ) SUSY QCD correction obtained 
by diagrammatic calculation 17 . 
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2. Two-loop SUSY QCD mass correction 

The pole mass nig of the gluino, which is de- 
fined by the complex pole s. p — (m§ — iTg/2) 2 of 
the gluino propagator, is given at the two-loop 
order as 



m g = M 3 + Sm g 



(i) 
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where the corrections Srn^ are expressed in 
terms of the one-loop and two-loop parts of the 
gluino self energy E(p) = T, K (p 2 )$ + T, M (p 2 ) as 



8m\ 

X (2 
0771- 



(1) _ 



-Re[M 3 E^ } (Ml 



-Re[Af 3 E^(M 2 ) + E$(M 2 )] 



+Re 



{M 3 E 



(i 



K (M 2 ) + ^(M 2 )} 

2v(l)/ 



x{E«(M 3 2 



2M£X K 



{Mi) 



+2M 3 E«(M 2 )} 



(2) 



The two-loop 0(a 2 ) correction Sm^ consists 
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8m^'' 2 \ where 



of two parts, 8m- = 8m 
8m- ' is the contribution of the diagrams with 

(2 2) 

only gluons and gluinos, while 8m y ~ ' is the re- 
maining contribution including quark and squark 
loops. Two-loop diagrams for these contributions 
are shown in Fig. ^ an d Fig- 121 respectively. In 
these figures, the wavy line, solid line without an 
arrow, solid line with an arrow, and dashed line 
with an arrow represent the gluon, gluino, quark, 
and squark, respectively. 
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Here M 3 is the running tree-level gluino mass. 
The dot in Eq. (0) denotes the derivative with 
respect to the external momentum squared p 2 . 

The SUSY QCD contribution to E(p) is gener- 
ated by loops with the gluino, gluon, quarks, and 
squarks. Masses and couplings in the lagrangian 
are renormalized in the DR scheme ^H] at the 
scale Q. Here we ignore SU(2)xU(l) breaking ef- 
fects in the loops, such as the quark masses and 
squark left-right mixings. This approximation is 
valid for the case where the gluino and squarks 
are sufficiently heavier than the quarks. Later 
we will briefly comment on the effects of these 
SU(2)xU(l) breakings. For simplicity, we also 
assume degenerate mass for squarks. 

The one-loop correction Sm^ in our approxi- 
mation is |15I16| 



8m- 



Cya s ( M 3 
= —J*, (5-6108-^. 

+—N q T F M 3 B 1 (Ml,0,mg), (3) 

where Cy = 3, Tp = 1/2, and N q — 6 is the num- 
ber of quarks. Parameters (a s , M 3 , m^) in Eq. 
© are the DR running ones at the renormal- 
ization scale Q. B\(p 2 ,mi,m,2) is the one-loop 
function in the convention of Ref. 20 . 



Figure 1. Two-loop 0(a 2 ) contributions to 

(2 1) 

8m- ' , without quark and squark propagators. 



(2 1) 

The contribution 5m~ ' is obtained by apply- 
ing the formula of the 0(a 2 ) QCD correction to 
the quark masses f2TJ in the DR scheme [22] to 
the SU(3) octet fermion. The result is 



(2,1) _ 



5m 



— j M 3 (-48 log— +36 log — 

+26 + 5tt 2 - 4tt 2 log2 + 6C 3 ) , 



(4) 



where ( 3 = Sn=i n 3 — 1-202. We have verified 
Eq. @ by explicit calculation of the diagrams. 
At Q = M 3 , the correction Q is 8mf' X) /M 3 ~ 
31(a s /7r) 2 - 0.03. 

(2 2). 

The contribution 8m- ' including quark and 
squark loops is calculated by decomposition into 
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Figure 2. Two- loop 0(a 2 s ) contributions to 

(2 2) 

Sm~ ' , with quark and squark propagators. 
Other diagrams obtained by charge conjugation 
are not shown. 



3. Numerical results 

We present some numerical results of the 0{a 2 s ) 
pole mass of the gluino, for the running tree-level 
mass M 3 (M 3 ) = 580 GeV which is close to the 
values in the SPSla point. The strong coupling 
constant within the standard model is given as 
a s (m z ) = 0.12. 

We first show, in Fig. the residual depen- 
dence of the one-loop pole mass and two- 

(2) 

loop pole mass m- on the renormalization scale 
Q, for the running squark mass m,q(Qo) — 800 
GeV at Qq — 580 GeV. All parameters in the 
formulas are evolved by 0(a 2 s ) renormalization 
group equations. For reference, the tree-level 
M 3 (Q) decreases from 589 GeV at Q = 400 GeV 
to 559 GeV at Q = 1400 GeV. We see that the Q 



two-loop basis integrals |23I24| using the integra- 
tion by parts technique |23I25I26| . and their nu- 
merically evaluation by the package TSIL |27) . 
We have analytically checked that the Q depen- 

(2) 

dence of 5m- is consistent with the two-loop 
renormalization group equation |14I28| of M 3 . 

(2 2) 

The explicit form of 5m~ ' is rather long. Here 
we just show, for reference, the form in the limit 
of vciq > M 3 : 

5mf 2 \m q » Af a ) = 
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The last term of Eq. (JSJ, which is independent 
of rriq, comes from the diagram (a) in Fig. [21 We 
have checked that the dependence of Eq. JSJ is 
consistent with the two-loop running of the gluino 
mass in the effective theory where squarks are 
integrated out 
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Figure 3. Dependence of the one-loop (dashed) 
and two-loop (solid) pole masses of the gluino on 
the renormalization scale Q. Mass parameters are 
(M 3 ,m 5 ) = (580,800) GeV at Q = 580 GeV. 



(2) 

dependence slightly improves by including Sm~ . 
One should however note that, contrary to naive 

(2) 

expectation, 5m- is much larger than the in- 
dependence of the one-loop result w^ 1 '- 
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Figure 4. The pole masses of the gluino at 
the one-loop with running masses (dashed), two- 
loop (solid), and modified one- loop mass with 
pole masses (dot-dashed), for the tree-level mass 
M 3 (M 3 ) = 580 GeV, as functions of m 4 (M 3 ). 



ings of squarks. Since these parameters break 
the SU(2)xU(l) gauge symmetry, their contri- 
butions to rrig should be suppressed by factors 
m^/m? or m 2 /m g compared to the gauge sym- 
metric contribution shown here. In addition, they 
only modifies contributions involving quarks and 
squarks in the third generation, while all gener- 
ations contribute to the SU(2) xU(l)-symmetric 

(2 2) 

part of dm- ' with equal weight. We therefore 
do not expect that these SU(2)xU(l)-breaking 
contributions are numerically relevant in future 
realistic studies of the SUSY particles. How- 
ever, detailed study in cases of light gluino and/or 
squarks is necessary for definite conclusion. 

4. Other two-loop corrections 

Beyond one-loop, there are also mass cor- 
rections involving Yukawa and electroweak cou- 
plings. For example, diagrams involving Yukawa 
couplings h q of the Higgs bosons/higgsinos to 
quarks and squarks in the third generation, such 
as those in Fig. give the 0(a s h 2 ) mass cor- 
rections. Analytic form of these contributions 



In Fig. 01 we compare m~ and m~ as func- 
tions of the running squark mass rriq(Q = M3). 

Here the renormalization scale is fixed at Q — 580 

(2) 

GeV. The two-loop correction 5m~ is positive 
and in the range of 8-15 GeV for m s = 400-1600 
GeV. This correction is 0(1 — 2) % of the one- 
loop result , with partial cancellation be- 
tween 5m~ 2,1 ^ and Sm^' 2 \ but still similar to or 
larger than the expected uncertainty in the mass 
determination at future colliders J2| • In Fig. ^ 
we also show the modified one-loop result where 
the running masses in Eq. are replaced by 
the corresponding pole masses. This modification 
corresponds to the inclusion of higher-order cor- 
rections by one-loop renormalization group equa- 
tions However, it is clearly seen that the 
resulting change of nig is much smaller than the 
two-loop correction Sm^K 

A recent work |30| writes down the explicit 
form of the 0(ct 2 s ) mass correction in terms of the 
basis integrals, including m q and left-right mix- 




Figure 5. Some diagrams of the 0(ct s h 2 ) contri- 
butions to the gluino pole mass. 



can be derived from general formulas of the two- 
loop corrections to the fermion pole masses j3U] 
in the approximation of massless vector bosons in 
the loops. As a numerical example, again in the 
SU(2)xU(l) symmetric approximation, the top 
Yukawa contribution for M3 = rriq = m^o = fi, 

tan/3 = 10 is 5mf M) /M 3 ~ (2, 10) x a s h 2 /(4Tr) 3 
for At = (M 3 , — M3), respectively, which is much 
smaller than the 0(a 2 ) contribution. We expect 
that the smallness of the Yukawa contribution 



Two-loop SUSY QCD correction to the gluino pole mass 



5 



would also hold in more general cases. 

5. Conclusions 

We have calculated the two-loop SUSY QCD 
contribution to the gluino pole mass, ignoring 
SU(2)xU(l) symmetry breakings in the loops. 
The 0(o? s ) correction to the gluino mass has been 
shown to be typically 1 — 2 %. For the case of 
M 3 (M 3 ) = 580 GeV, this correction is similar to, 
or larger than, the expected uncertainty in the 
mass determination from precision measurements 
at future colliders. The two-loop correction would 
be therefore important in the extraction of M3 
from experimental data and the determination of 
the SUSY breaking at the unification scale. 
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